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Abstract 

A analysis of several SUSY GUTs recently discussed in the literature 
is presented. We obtain global fits to electroweak data, which include gauge 
couplings, gauge boson masses, BR{b sj) and masses of fermions of all 
three generations and their mixing angles. Thus we are able to test gauge 
unification, radiative electroweak symmetry breaking, SUSY sector ( - in the 
context of supergravity induced SUSY breaking) and the Yukawa sector in each 
particular model self-consistently. One of the models studied provides a very 
good fit with ^ ^ for 3 degrees of freedom, in a large region of the allowed 
SUSY parameter space. The Yukawa sector works so well in this case that the 
analysis ends up testing the MSSM constrained by unification. Adopting this 
point of view, in the second part of this talk we focus on the details of the fit 
for BR{b 57) and discuss the correlations among Sm^^^^ , as{Mz) and a 
GUT threshold to as{Mc)- We conclude that an attractive S'O(10)-derived 
regime of the MSSM remains a viable option. 
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A analysis of several SUSY GUTs recently discussed in the literature is pre- 
sented. We obtain global fits to electroweak data, which include gauge couplings, 
gauge boson masses, BR{b S'y) and masses of fermions of all three generations 
and their mixing angles. Thus we are able to test gauge unification, radiative elec- 
troweak symmetry breaking, SUSY sector { - in the context of supergravity induced 
SUSY breaking) and the Yukawa sector in each particular model self-consistently. 
One of the models studied provides a very good fit with ~ 1 foJ" 3 degrees of 
freedom, in a large region of the allowed SUSY parameter space. The Yukawa 
sector works so well in this case that the analysis ends up testing the MSSM con- 
strained by unification. Adopting this point of view, in the second part of this 
talk we focus on the details of the fit for BR{b — » 57) and discuss the correlations 
among 5m^^^^ , as{Mz) and a GUT threshold to oia{MQ). We conclude that an 
attractive 5'O(10)-dcrived regime of the MSSM remains a viable option. 

1 Introduction 

Reaching beyond the Minimal Supersymmetric Standard Model (MSSM) we 
apply global analysis to the search for grand unified theories (GUTs). That 
enables us to compare specific GUT models which become the MSSM as effec- 
tive theory below the unification scale. In practice, it means that we test gauge 
unification, radiative electroweak symmetry breaking and the SUSY sector in 
the same way as the MSSM analysis constrained by unification EJ, but in addi- 
tion we also test the Yukawa sector of the theory versus the observed fermion 
masses of all three generations and the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. Clearly, we agree with the MSSM analysis for models which 
describe the Yukawa sector very well. 

There are plenty of models to be tested in this kind of analysis. For the 
beginning we start with models based on SO(IO) gauge symmetry. SO(IO) 
SUSY GUTs have been recognized as excellent candidates for an effective field 
theory below the Planck/string scale. They maintain the successful predic- 



1 



tion for gauge coupling unification and provide a powerful and very economic 
framework for theories of fermion masses. This is because all the fermions of 
one generation are contained in the 16 dimensional representation of SO (10) - 
thus fermion mass matrices are related by symmetry. In the most predictive 
theories, the ratio of Higgs vevs - tan/3 - is large, and the top quark is naturally 
heavy as found experimentally. However, as a consequence of large tan/3 there 
are potentially large supersymmetric one-loop effects at the weak scale which 
could play an important role in fitting the fermion masses and mixings, and the 
FCNC processes like the observed 6 — > 57 decay rate. Thus a self-consistent 
analysis becomes more powerful (and restrictive) than in a low tan/3 scenario. 

In sections 3 and 4, we present the results of such a complete top-down 
analysis. In section 5, we focus on the MSSM constrained by the best working 
GUT model and analyze it similarly to generic MSSM analyses constrained by 
unification. ^ We present the best fit correlations among various contributions 
to the 5 — > 57 decay rate, SUSY corrections to m?,, and rather low values of 
as{Mz) - all in the (toq, M1/2) SUSY parameter space. The phenomenological 
implications are discussed as well. 

2 Global Analysis 

Details of our numerical procedure are described ini. The analysis starts at 
the GUT scale AIq, which is a free parameter itself, with unified gauge cou- 
pling ac, Tiy free parameters entering the Yukawa matrices^, and with £3 as 
a one loop GUT threshold correction to as{Mc) |. We assume supergravity 
induced SUSY breaking. At the scale Mg we introduce standard universal soft 
SUSY breaking parameters mo, M1/2 and Aq, and non- universal Higgs masses 
mn^ and mn^ ■ The fi parameter and its SUSY breaking bilinear partner B 
are introduced at the Z scale, since they are renormalized multiplicatively and 
do not enter the RGBs of the other parameters. At the Z scale we match the 
MSSlVLdirectly to SU{S)c x U{l)em, thus leaving out the SM as an effective 
theoryB. Electroweak symmetry breaking is established at one loop in the pro- 
cess of the minimization by fitting the observables in table |l|. Within the 
MSSM, we calculate one- loop corrected W and Z masses and G^, corrections 

" However, our MSSM analysis is simplified. We give up on most of the precision electroweak 
observables and keep only ten observables of the MSSM analysis (see the next section for 
details). We believe that this reduction does not bias our results in any significant way since 
the asymmetries and Z lineshape parameters do not present dominant constraints, especially 
if most of the SUSY particles are rather heavy (as happens for the best working model of 
our analysis). 
Clearly, riy is model dependent. 

Mq is defined as the scale where the gauge couplings ai and 02 are exactly equal within the 
one-loop GUT threshold corrections. By ac we actually mean the value ai{MG) = 0:2 (Mq) 
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Table 1: Experimental observables of the global analysis. 



Observable 


Central 


(T 


Observable 


Central 








value 








value 




1. 


Mz 


91.186 


0.46 


11. 


Mb - Me 


3.4 


0.2 


2. 


Mw 


80.356 


0.40 


12. 




180 


50 


3. 




1.166 • 10-5 


1.2 • 10-^ 


13. 


md/ms 


0.05 


0.015 


4. 




137.04 


0.69 


14. 


0.00203 


0.00020 


5. 


as{Mz) 


0.118 


0.005 


15. 




105.66 


0.53 


6. 


Mt 


175.0 


6.0 


16. 


Me 


0.5110 


0.0026 


7. 


mh{Mb) 


4.26 


0.11 


17. 


Vus 


0.2205 


0.0026 


8. 


Mr 


1.777 


0.0089 


18. 


Vcb 


0.0392 


0.003 


9. 


Pnew 


-0.6 • 10-^ 


2.6 • 10-^ 


19. 


VublVcb 


0.08 


0.02 


10. 


B{b sj) 


2.32 • 10-"* 


0.92 • 10-* 


20. 


Bk 


0.8 


0.1 



to the p parameter from new physics outside the SM, and the amplitude for the 
process 6 — > 57. When crossing the Z scale, we compute the complete one loop 
threshold corrections to as and a, whereas only those one loop threshold cor- 
rections to the fermion masses and mixings enhanced by tan/3 are computeds. 
The .Z-scale amplitude for 6 ^ 57 is matched to the coefficient C'j{Mz) which 
is then renormalized down to the scale Mb, 

o 8 

C'jff{Mb) = r^"^^C,{Mz) + ^ (^^ - ryM) C^^Mz) + C2{Mz) Kv""' , (1) 

1=1 

based on mixing of the electromagnetic operator with the chromomagnetic and 
current-current operators (coefficients Cs and C2) in the leading log approxi- 
mation; with 1] = as{Mz) / (y-s{Mb). Finally, the branching ratio 

BR{b sj) = ^-^^^^ \C'r"i^^b)\' BR{b ^ cev) , (2) 

\VcbY TTg[Mc/Mb) 

where BR{b cev) = 0.104, a ~ 1/132.5, the phase-space function giz) — 
1 — 8z^ -(-8z^ ~ ~ 24z* log 2, and hiS and a^'s in eq.(|l|) are given in refB. The 
values of the CKM matrix elements and quark masses in (0) are consistently 
calculated in the actual fit within a particular model. 

Note that Mz, Mw, a, and the lepton masses are known so well that 
we have to assign a theoretical error as their standard deviation ( — the cor- 
responding cr's are underlined in table [|). We estimated conservatively the 
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theoretical error to be 0.5% based on the uncertainties from higher order per- 
turbation theory and from the performance of our numerical analysis. The 
error on is estimated to be within 1% due to the fact that in addition to 
the uncertainties mentioned above we neglect SUSY vertex and box corrections 
to Ar. Also note that e^c, the observable of CP violation, has been replaced by 
a less precisely known hadronic matrix element Bk ■ Similarly, the light quark 
masses are replaced by their ratios (at the scale IGeV), and mc{Mc) by the 
difference Mj, — Mc, since the latter quantities are known to better accuracy. 

= (m^ — mfj/rri^ is the Kaplan-Manohar-Leutwyler ellipse parameter. 
Finally note, that the data in table ^ are divided into two groups. The first 
ten, i.e. five observables in the gauge sector {Mz, Mw, G^^ a and as), masses 
of the third generation fermions, Pnew and BR{b — > S7) can also be included 
in tests of the MSSM constrained by unification with no particular underlying 
GUT. The other ten correspond to six light fermion masses and four indepen- 
dent parameters of the CKM matrix and test primarily the Yukawa sector of 
a GUT model. 

In addition, the function is increased significantly by a special penalty 
whenever a sparticle mass goes below its experimental limit. 

3 Results for SO(IO) GUT Models with Four Effective Operators 

The analysis, as described above, was used to test simple SO(IO) models. First 
we checked the performance of the nine models with four effective operators 
suggested incl. The models were defined by a unique choice of the operators 
O33 = I63 10 I63 and O12 = I61 10 {^f I62 at the GUT scale (A's stand 
for adjoint states and S for singlets). There were six 022 operators which all 
gave the same 0:1:3 Clebsch relation between the 22 elements of the Yukawa 
matrices for the up quarks, down quarks and charged leptons, thus introducing 
Georgi-Jarlskog relation in all models. Finally, the models were disti|nguished 
by a choice of the O23 operator. Nine operators were suggested in refQ leading 
to different predictions at low energies. 

Our analysis has shown that the choice of the 023 operator is indeed 
significant (see fig. la) and that model 4 is by far the best working model. In 
this case, O23 = I62 (^) 10 (^) I63. Note that with four effective operators we 
have Hy = 5 free parameters in the Yukawa matrices at the GUT scale. That 
leaves the function (out of the 20 observables in table |l|) with 5 degrees of 
freedom (d.o.f.). We show in fig. lb that the performance of model 4 does not 
get significantly better in a larger SUSY parameter space. We checked that 
the same is true for the other models. We also checked that no substantial 
improvement of the performance of model 4 can be achieved by neglecting one 
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Figure 1: Results of global analysis for SO(IO) GUT models: (a) nine models of section 3; (b) 
model 4 [solid (dashed, short-long-dashed, dotted) lines represent contour lines of constant 
= 16 (15, 14, 13) per 5 d.o.f.]; (c) model 4c [solid (dashed, dotted) lines represent contour 
lines of constant = & (3, 1) per 3 d.o.f.]. 



out of the twenty observables given in table QEI. On the other hand, we have 
found that a significant improvement is possible by adding one new operator, 
contributing to the 13 and 31 elements of the Yukawa matrices. 

4 Results for SO(IO) GUT Models with Five Effective Operators 

Next, we analyzed two models recently derived from the complete SO(IO) 
SUSY GUTsEl. The models were constructed as simple extensions of model 
4. Different label (a,b,...f) refers to the different possible 22 operators. In 
the extension to a complete GUT the different 22 operators lead to inequiv- 
alent theories due to different U(l) charge assignments. When one demands 
"naturalness" , i.e. includes all terms in the superpotential consistent with the 
symmetries of the theory one finds one and only one 13 operator for models 
4a and 4c. The 22 and 13 operators of model 4c are O22 = I62 I 10 ^ I62 

and ©13 = I61 {^Y 10 ^ I63 • With the 13 operator ny—7, which implies 
3 degrees of freedom. The results of the global analysis are given in fig.lc. 
Similarly to fig. lb, these figures show the contour plots of the minimum 
in the {mo, AI1/2) plane, for three different fixed values of the parameter 
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fi{Mz) — 80, 160, 240 GeV. All initial parameters other than {mo, M1/2, n} 
were subject to minimization. 

The fits get worse as fj, increases because the SUSY corrections to fermion 
masses and mixings increase with fj,. As n gets larger they can only be kept 
under control by larger squark masses. Varying fj, freely actually results in its 
approaching the lowest possible value. The lower bound on fi{Mz) is deter- 
mined by the chargino mass limit from direct searches and is correlated with 
A/1/2 ■ When the value of /x is fixed, as in fig's Ib-c, the chargino mass limit 
then sets a sharp lower bound on -M1/2, which is explicitly visible in each of 
the figures. Figures lb and Ic were obtained assuming m^- > 65GcV. In 
our further analysis ii{Mz) has been fixed to llOGeV and m^- > 85GeV was 
imposed. Figures 2a-b show explicitly that the structure observed in fig.lc 
originates from the two distinct fits corresponding to two separate minima of 
the global analysis. The fits are primarily distinguished by the sign of the 
& ^ 57 decay amplitude. ( See section 5 for more details.) We do not show 
the results for negative values of /i. In this case, the chargino contribution to 
6 — > S7 interferes constructively with the already large enough SM and charged 
Higgs contributions. As a result, the fits get much worse, with well above 
10 per 3 d.o.f., and that disfavors this region of the SUSY parameter space. 
Similar observations were also made by W. HoUik et alu. 

Model 4a is defined by different 22 and 13 operators and gives the fits with 
the best —4-6 in most of the parameter space. (It yields also ~3, hx± 
only for the corner in the SUSY parameter space with large mo, and /x.)0 

Whether or not these particular models are close to the path Nature has 
chosen remains to be seen. One important test will be via the CP violating 
decays of the B meson. Models 4c and 4a both predict a narrowly siaread value 
sin2a ~ 0.95, whereas in the SM the value of sin2Qj-|is unrestrictedB. Another 
important test may come from nucleon decay rates cl. 



5 MSSM Analysis Constrained by the Best Working GUT Model 

Since the fermionic sector of model 4c works very well, we can regard the 
analysis in this case as a test of the MSSM with large tan/3. Yet, some model 
dependence is still present. It is, first of all, the introduction of £3, a GUT 
threshold to a^. It is actually the only GUT threshold introduced in our 



study. Note in subsection 5^ that non-zero £3 is imposed by the low energy 
phenomenology rather than by physics at the GUT scale. 

For the Yukawa matrices, the exact equality of the 33 elements is assumed. 
The remaining Yukawa entries are small and decouple from the MSSM RGEs 
for the gauge and third generation Yukawa couplings and diagonal SUSY mass 
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parameters. Thus they have no effect on the calculation of the Z-scale values 
for the first nine observables in table |l|. The BR{b — > 57) is affected by some of 
these entries. Model dependence comes from the chargino contribution which 
contains inter-generational cl-Il squark mixing. That contribution is tan/3 
enhanced which makes it non-negligible. The mixing is completely induced by 
the off-diagonal entries of the Yukawa matrices in the RG evolution. Next, 
note that the light quark and lepton masses and CKM elements do not exert 
any significant pull on the best x^'s of model 4c □. We conclude that our 
results presented in this section are not sensitive to the structure of the Yukawa 
matrices except for the model dependent 23 mixing which is significant for the 
BR[b ^ sj). 

5.1 Results for BR{b 57) 

We find that the b ^ sj amplitude is dominated by the SM, charged Higgs, 
and tan/3 enhanced chargino contributions. To understand our results, let's 
estimate first what to expect from the SUSY contribution to C7{Mz) and define 
the Z scale ratios for Higgs and chargino contributions = Cy^^Cy'^*^'' 
and rC^) = c\^^ /c\^^''\ Equation (0) then reads 

^(MSSM) eff ^ ^la ^(SM) ^ ^(«) ^ ^(C)) ^ ^ ^^^a. (3) 

1=1 

where we used C2{Mz) = 1 and Cs "C Cy, C2. Since it is well known that 
(j{SM) eff (-j(SM) _|_ ^8^^ hif]"-^) would yield about the right value for 

the BR{b 57) we infer that either 

^(C) ^ _^iH) ^ ^(M5SM) eff ^ ^ ^(5M) eff 

or 

r(C) „ _^iH) _ 4 60 , for C^/^^SM) eff ^ _ ^(SM) eff ^ 

For the last estimate, the numerical results C^'^*''^' = —0.190, = 0.679, and 
Si=i hi'q"'^ = —0.168 were used — computed for as{Mz) ~ 0.118. 

The charged Higgs contribution always interferes constructively with the 
SM contribution. We get < r'^' < 1.3, depending on the mass of the H~ . In 
the first case, especially if A^'' and r*^*-^^ are non-negligible, eq.(|^) means that 
the chargino part must interfere destructively^ with the SM and charged Higgs 



"^Since cq's ^ and are valid only approximately, the case 1 + -l-rC-^) ?a 1 in principle 
also allows for a constructive interference < r^^\ r^H) ^ 1 ti^g region in parameter 
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contributions, practically cancelling the latter. The enhancement by tan/3 of 
the chargino contribution has to be compensated for by rather large sparticle 
masses. In the second case, described by eq.(||), large destructive chargino 
interference is required to outweigh the combined SM and H~ contributions 
and to flip over the overall sign of the amplitude. Quite amazingly, it is not 
so difficult to arrange since the chargino contribution is the only one enhanced 
by large tan/3. However, large sparticle masses obviously suppress the effect. 
The lesson is that we can expect the two cases to work in a complementary 

SUSY parameter space and have sign (^j^^'^'^*^) =(-) sign cf'^^ for the best 
fits in the region with large (low) mp and/or M1/2, respectively. 

These expectations are indeed realized in the best fits of model 4c. Figure 
2a (2b) corresponds to the case when the 6 ^ 57 amplitude in MSSM is of 
the same (opposite) sign as the SM amplitude. As anticipated, these two cases 
work in complementary regions of parameter space. Figures 2c-d show how 
well the fits describe the measured value of the BR{h 57). In these figures 
(and similarly in the following ones) we show the contour lines of fig's 2a-b 
in the background for reference. As can be seen, the BR{b — *■ 57) presents a 
major constraint since whenever the values go up, the agreement with the 
observed 6 — > 57 decay rate gets worse. For the flipped sign of the amplitude 
it is interesting that the lightest stop ti can be sufficiently light (see fig.2f) to 
maintain the required large size of the chargino contribution even for mo at 
2TeV. That can be arranged due to significant left-right stop mixing and the 
freedom in varying the trilinear scalar coupling Aq{Mc)- The corresponding 
tL-tR mixing angle changes between 31° to 40° across the (mo,Mi/2) plane 
while the trilinear coupling At{Mz) ranges from —300 to — ITOOGeV in this 
case. Figure 3 shows the contour plots of constant r^^' and for different 
signs of C'^^'SSAf^ Qj-^g compare the numerical results in figures 3a vs. 3c 
and 3b vs. 3d with the approximate relations (^) and (|^). 

There is one striking feature which is common to both cases. It is that 
both fits would like to have the BR{bj-^ 57) below rather than above the 
CLEO experimental value 2.32 x lO""*. B 

In the first case, the tan/3 enhanced chargino contribution tends to be too 
large when going against the charged Higgs contribution, since the latter is not 
tan/3 enhanced. The fit clearly tries to make the Higgs contribution as large as 
possible (see fig. 3a). As a result, the charged Higgs (and then also the whole 
Higgs sector) tends to be light. We get, for instance, the best fit value of the 
pseudoscalar mass uia < lOOGeV everywhere in the (mo,Mi/2) plane. 

space where rrijj- and sparticle masses are large. This option, however, does not result from 
our best fits, as already mentioned in the discussion on negative /j parameter in the previous 
section. 
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Figure 2: Model 4c global analysis results for the two possible signs of Cf''^'^'^ . The best 
fit contour plots of (a),(b) [solid (dashed, dotted) lines correspond to = 6 (3, 1) per 3 
d.o.f.]; (c),(d) BR(6 — » s-y) X 10* ; (e),(f) the Ughtest stop mass m^^ [in GeV]; with the x^ 
contour lines in the background. 
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Figure 3: Relative contributions to C7 for the two possible signs of C^^^^^' . Contour plots 
of (a),(b) r'-'^'; (c),(d) r^'-^^; {e),(f) j-(C'23). ^j^g ^2 cqji^quj- lines in the background. 
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with the contour lines in the background. 
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In the second case, when the sign of C7 is flipped by the chargino contribu- 
tion, this contribution tends to be not big enough, especially if > 300GeV. 
Hence & — > 57 would prefer a heavy charged Higgs. However, X)ne cannot have 
good fits with the Higgs sector much heavier than squarkso and thus 
varies quite a bit across the SUSY parameter space (see fig. 3b). When M1/2 
gets below SOOGeV, 6 — > 37 is no longer a strong constraint and two separate 
minima can be found in the course of optimization. The two fits work equally 
well (x^ in each case stays below 1 per 3 d.o.f.) and differ only by the best fit 
values of the Higgs masses. One minimum corresponds to m^o and vtla settled 
at the experimental lower limit (65GeV, in our analysis) while these masses 
gradually rise up to 700GeV in the second "valley" . In the allowed corner with 
mo < 700GeV, the two "valleys" approach each other and eventually coincide. 
The effect of the doubled minimum is indicated in figures 3b, 3d, 3f and 4b 
with the solid black (dashed gray) contour lines corresponding to the heavier 
(lighter) Higgs sector. 

Can the NLO QCD corrections to BR(6 57) improve the fits? The 
calculation of these-corrections has not been completed in the MSSM. Here 
we just summarizill that the rate goes up if the unknown MSSM match- 
ing corrections satisfy c^^^'^'^'''^^'' j c^^'^'^''''-^'' ^ —1, (fc = 7,8), in the case 
when sign q^^'^^^^'^^ — _|_ gigp, . Thus the allowed SUSY parameter 

space is likely to increase in this case. In the second case, with the oppo- 
site sign of the 6 — > 37 amplitude, the NLO correction will more likely make 
the fit worse. In this case, the net NLO correction increases the rate only if 

IjiMSSM) (1) l(j(SM) (1) < _7 S 

In summary, if the future experimental analysis confirms the_discrepancy 
between the CLEO measured value and the NLO SM calculationB, the MSSM 
with large tan/3 could be the solution. Similarly, if the NLO calculation is com- 
pleted for the MSSM, and if it turns out to enhance the LO result as occurred 
for the SM, then the large tan/3 regime will apparently have no problem fitting 
the 6 — > 57 rate exactly. 



5.2 Role of c-t mixing for BR{b 57) 

In analogy to previously defined r^^^ and r^^'^ we introduce r^*^^^^ = cj*^^'^'' / C^" 
where C) is the CL-th mixing contribution to the coefficient C7 at the scale 
Mz- We show the contour plots of the constant r'^'-^^^-' in fig's 3e-f. While 
the dominant chargino contribution is that proportional to the tL-tn mixing, 
^jC23) ]-,gj.Qj^gg iniportant because of the destructive character of the interfer- 
ence among the partial amplitudes. As one can see, the c^-i^ mixing term can 
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be comparable with the SM and H contributions. More importantly, it al- 
ways interferes constructively with them. In the case when sign (^(^^'^■S'^^) _ ^ 
sign Cy'^*^'' , this term helps to counterbalance the large contribution of the left- 
right stop mixing. In the complementary case, when the chargino contribution 
flips the sign of the net amplitude the CL-tL mixing makes it more difficult to 
happen. As a result, it has different consequences for the two fits in figures 2a 
and 2b, especially important in the region (mo < lOOOGeV, Afi/2 ~ 350GeV) 
where the fits start getting worse. These observations are, however, model 
dependent and their validity relies on the boundary conditions assumed at the 
GUT scale. 



5.3 Correlations among £3, as{Mz) and 6m^^^^ 

The results on 6 — > 37 show that it is important to have a destructive interfer- 
ence among the partial contributions to & — > 57. However, with the universal 
boundary conditions at the GUT scale that can be arranged only for a specific 
sign of the fi parameter: /i > in our conventions. That in turn correlates with 
the positive sign of the SUSY corrections to the b quark mass. This fact implies 
strong constraints on the SUSY parameter space from nih since 5mf^^^ gets 
a dominant contribution from the tan/3 enhanced gluino exchange. The giuino 
correction can be explicitely suppressed by heavy squarks and by the chargino 
correction which enters with the opposite sign. It can also be reduced by lower 
values of as{Mz) because of the strong coupling in the vertices of the gluino - 
6-squark diagram. The contour plots of the net SUSY correction to mi,{Mz) 
in the best fits are shown in figures 4a-b where one can see that these effects 
are quite effective in reducing 5mf^^^ . 

Note that the effect of large positive 5m^^^^ is also reduced by a lower 
value of as{Mz) indirectly — due to the RG evolution of the current b mass 
from Mz down to Affc. In our analysis, which assumes larger uncertainty for 
as than for mi,{Mi,), it pushes as{Mz) down (see figures 4c-d). One can trade 
lower values of as for higher values of m?, provided a smaller uncertainty a (as) 
is assumed. Rather low values of as{Mz) are, however, difficult to obtain from 
an exact gauge coupling unification. Our analysis shows that a few per cent 
negative correction to as{Mo) is enough to get as{Mz) ^ 0.118. The best fit 
values of €3 are presented in figures 4e-f. Clearly, as squarks get lighter the 
SUSY correction to mb{Mz) (fig's 4a-b) has to be increasingly reduced by the 
lower values of as{Mz) (fig's 4c-d), which in turn requires a more substantial 
departure from the gauge coupling unification (fig's 4e-f). 63 can be genemted 
by the spread in masses of heavy states integrated out at the GUT scale Q. 
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6 Conclusions 



At the present time, when direct evidence for physics beyond the SM evades ex- 
perimental observations, global analysis serves as the best test of new physics. 
Our results show that minimal SUSY SO(IO) models remain among the candi- 
dates for the theory below the Planck scale and that they can provide specific 
guidelines in the search for signatures of new physics at the electroweak scale. 
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